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ABSTRACT 

Context. For many years, there seemed to be significant differences between the continuum intensity distributions derived from 
observations and simulations of the solar photosphere. 

Aims. In order to settle the discussion on these apparent discrepancies, we present a detailed comparison between simulations and 
seeing-free observations that takes into account the crucial influence of instrumental image degradation. 

Methods. We use a set of images of quiet Sun granulation taken in the blue, green and red continuum bands of the Broadband 
Filter Imager of the Solar Optical Telescope (SOT) onboard Hinode. The images are deconvolved with Point Spread Functions (PSF) 
that account for non-ideal contributions due to instrumental stray-light and imperfections. In addition, synthetic intensity images are 
degraded with the corresponding PSFs. The results are compared with respect to spatial power spectra, intensity histograms, and the 
centre-to-limb variation of the intensity contrast. 

Results. The intensity distribution of SOT granulation images is broadest for the blue continuum at disc-centre and narrows towards 
the limb and for longer wavelengths. The distributions are relatively symmetric close to the limb but exhibit a growing asymmetry 
towards disc-centre. The intensity contrast, which is connected to the width of the distribution, is found to be (12.8 ± 0.5) %, (8.3 ± 
0.4) %, and (6.2 ± 0.2) % at disc-centre for blue, green, and red continuum, respectively. Removing the influence of the PSF unveils 
much broader intensity distributions with a secondary component that is otherwise only visible as an asymmetry between the darker 
and brighter than average part of the distribution. The contrast values increase to (26.7 ± 1.3) %, (19.4 ± 1.4) %, and (16.6 ± 0.7) % for 
blue, green, and red continuum, respectively. The power spectral density of the images exhibits a pronounced peak at spatial scales 
characteristic for the granulation pattern and a steep decrease towards smaller scales. The observational findings like the absolute 
values and centre-to-limb variation of the intensity contrast, intensity histograms, and power spectral density are well matched with 
corresponding synthetic observables from three-dimensional radiation (magneto- )hydrodynamic simulations. 

Conclusions. We conclude that the intensity contrast of the solar continuum intensity is higher than usually derived from ground- 
based observations and is well reproduced by modern radiation (magneto-)hydrodynamic models. Properly accounting for image 
degradation effects is of crucial importance for comparisons between observations and numerical models. 
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■ 1. Introduction granular intensity fl uctuation (see, e.g. jDeubner & Mattigll975l : 

T~ '„ , . , . , . . . , , iMuller et al1 l2007). Its value is straightforward to calculate and 

> ; Modern numerical simulations of convection in the upper solar ft ^ d j connected t0 the basic mec hanisms of convection as 

. - , atmosphere achieve a high degree of re a^For example in h . g sensitiye tQ ±e distribution of bri ht (hot) and dark (cool) 

?S • land T Paper ' I Stein&Nordlund! <|200rJ show that their surface stmctures It is there fore troubling that contrast values 

H ] MHD simulations are very successful at reproducing a number deriyed from simulations do not agree with observed contrast 

of observed properties of granulation. Other remarkable exam- yalues e .g.. iNordlundl 1984: Keller 2006) with synthetic 

pies are tlie exc d ent agreement betw een observed and synthetic cQntrast yaluefj being usually significant i y higher . For example, 

Fe line profiles jAsplund et alj | 2000|) and between i objgryedand | Uitenbroek et ^ J2007I) compared G-band images of the quiet 

pl ated magnetic faculae dCarlsson et al. | |20Q4j; i KelleretalJ Sun Th contrastfj of R j % for ^ (reconstmcted) ob _ 

|2QQ4i | Steiner| | 2005[) . As these diagnostic tests demonstrate the seryed image and 2L5 % for ^ Aedc { from a model 

fidelity of the simulations, they can be used to achieve a deeper , Stein & Nordlund ,2006). 
understanding of the physical mechanisms acting in the solar at- 
mosphere and provide access to the un-observable subsurface 

regions. In recent years, extracting observables from Simula- Qualitatively, empirical studies show that the granulation 

tions by synthesized spectra and intensity maps has grown into a contrast decreases with increasing wavelength and decreases 

powerful tool to u nders tand and interpret observations (see, e.g., with_ increasing distance from the disc-center. I Wilken et al. | 

Khomenk^itiDllQOl; lLangangen et alll2007l and many more). <H297D used the German Vacuum Tower Telescope (VTT) to ob- 

Despite their success, there exist a number of diagnostic tests tam granulation contrast at different hehocentac positions at a 

that show considerable discrepancy. One notorious example is wavelength of 550 nm. They derive (13.5 + 1) % for fi = 1.0 

the (rms) contrast of the continuum intensity emerging from the (disc-centre), 8-9 % at /i = 0.1, and down to 3.8 % very close 

low solar photosphere, also referred to as granulation contrast or to the solar llmb at A* = 0-07. The wavelength dependence and 

center-to-limb variation of the contrast provide a stringent di- 

Send offprint requests to: sven.wedemeyer-bohm@astro.uio.no agnostic test for simulations as it connects to the temperature 

* Marie Curie Intra-European Fellow of the European Commission stratification of the photosphere. 
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Fig. 1. SOT/BFI observations at disc-centre: Close-ups of intensity images in the blue (left), green (middle), and red (right) contin- 
uum band. 



In observations many factors affect the contrast and most 
of them are poorly understood. Imperfections in the optics and 
straylight significantly reduce the contrast but are very difficult 
to measure and correct for. In addition, image degradation by the 
Earth atmosphere (seeing) has a huge impact on the image qual- 
ity and can at best only partly be corrected for. The determination 
of the centre-to-limb variation of the mean abso lute continuum 
intens ity already requires enormous efforts (e.g., Nec kel & Labsl 
1984). Owing to these hardly controllable factors, a large range 
of contrast values are reported in the lit erature. See for ex- 
ample the review by [S pruit et al .l (1 19901) and mor e recently 
Table 2 by San chez Cuberes et al.l (120001) and Fig. 1 bv lKiselmanl 
(2008). Continuous advancement of observational methods, in 
particular the development of adaptive optics and image re- 
construction techniques, lead to an increasing trend of the re- 
porte d contrast values. These methods like for example spe ckle 
(e.g. Jvon der Luehell993l:lde Boerll993l:IW5ger et alj2008l) and 
Phase Diversitv/MOMFBD (Ivan Noort et al.l2005l) involve mea- 
surements of the effective point-spread-function (PSF). Despite 
the high level of sophistication and impressive results, a com- 
plete reconstruction (i.e., full removal of the combined effect 
of seeing-induced degradation and instrumental aberrations) can 
never be obtained. Contributions to the PSF that are important 
for accurate photometry - essential for the determination of the 
granulation contrast - are often hardly known. This is in par- 
ticular true fo r stray-light produced inside the optical instrument 
(Matt ig|1983l) . Already a simple uniform intensity offset reduces 
the resulting intensity contrast. In most cases, however, the stray- 
light contribution is likely to be anisotropic and varying over the 
field of view (FOV). 

With the advent of the Solar Optical T elescope (SOT 



Tsuneta et al.l 120081; 



Ichimoto et al.l 120081; ISuematsu et al.l 
2008; Shimiz u et al.l 120081) onboard the Hinode spacecraft 
dKosugi et al.l 12007). the solar community has now access to 
space observations that have sufficient spatial resolution to 
resolve the scales that are important for an accurate determina- 
tion of the granulation contrast. A space-born instrument has 
the obvious advantage of the absence of seeing so that only 
the instrumental PSF or stray-light needs to be considered. 
IWedem ever-Bohml J2008L hereafter referred to as Paper I) 
derived a set of PSFs for the combined effect of the Broadband 
Filter Imager (BFI) and the Optical Telescope Assembly (OTA) 
of SOT. He concluded that there is a small but significant 
amount of stray-light, although the instrument still operates 



close to the diffraction limit. The effect on the intensity dis- 
tribution, including the granulation c ontrast, could t hus b e 
significant. This view is supported by iDanilovic et al.l (|2008), 
who find that the difference in contrast between observations 
with the Spectro-polarimeter (SP) of SOT and state-of-the-art 
numerical simulations is strongly reduced, when taking into 
account instrumental image degradation. 

In this paper, we exploit the extensive Hinode data base that 
has been built up through its successful operations. Over time, 
the various observing programs have covered the quiet Sun at a 
whole range of observing angles providing the data to facilitate a 
robust statistical study of the wavelength dependence and center- 
to-limb variation of granulation contrast. Furthermore, we con- 
sider different independent numerical models. The paper is or- 
ganised as follows: the observations are described in Sect. 12.11 
followed by a description of the PSFs in Sect. l2.2l and the method 
of deconvolution in Sect. 12.31 The numerical simulations and the 
intensity synthesis are introduced in Sect. [3] The observations 
and the simulations are compared in detail for disc-centre im- 
ages in Sect. [4] followed by an analysis of the centre-to-limb 
variation of the intensity distribution and the continuum intensity 
contrast in Sect. [5] Finally, the observations and simulations are 
compared in terms of spatial power spectral density in Sect. [6] 
Discussion and conclusions can be found in Sects. [7] and [8] re- 
spectively. 



2. Observations 

2.1. SOT filtergrams 

The observations analysed in this study were obtained with 
the Broadband Filter Instrument (BFI) of the Solar Optical 
Telescope (SOT) onboard the Hinode satellite (see Sect. 1 for 
references) in the period from November 2006 until February 
2008. Wide -band filtergrams in the blue, green and red contin- 
uum from various dates are considered (see Fig.[TJ. The nominal 
central wavelengths of these filters are 450.45 nm, 555.05 nm, 
and 668.40 nm, respectively. The transmiss ion profiles of thes e 
filters, which are available in Solar Soft dFreeland et al.l l2000). 
have FWHMs of 0.22 nm, 0.27 nm, and 0.31 nm, respectively. 
As in Paper I, the position of the FOV and individual subre- 
gions on the solar disc is here reduced to the heliocentric po- 
sition fi - cos 8 with 8 being the observing angle. It is justified 
by the finding that the intensity profiles across the solar limb 
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do not show noticeable differ ences between different pos itions 
along the limb (N/W/S/E) (cf. lLanghans & Schmiddl2002l) . The 
FOV is 223 " x 1 12 "for images that use the whole detector and 
112" x 112 "for images that use only half. The pixel scale of the 
detector of only 0'.'054 is so much smaller than t he spatial reso- 
lution of the optical system (0'.'2 - 0'.'3, see, e.g. JShimizu et al.1 
2008) that images with 2x2 pixel-binning can be used next to 
unbinned images without causing significant differences in terms 
of contrast and intensity distribution. A total number of 584 im- 
ages of quiet Sun granulation are analysed, covering heliocentric 
positions from disc-centre to the solar limb. The images were se- 
lected to contain quiet Sun only, active region observations were 
excluded. Basic data processing (i.e., dark subtraction and flat- 
fielding) was preformed using routines from the Hinode branch 
of Solar Soft. 



2.2. Point Spread Functions 

We use the point spread functions for the Solar Optical 
Telescope as they are calculated in Paper I. There are individ- 
ual ideal and non-ideal PSFs for the blue, green, and the red 
continuum band. The ideal PSFs represent the image degrada- 
tion due to diffraction at the main aperture stop of the tele- 
scope, which has a diameter of 50 cm. The central obstruction 
and the spiders are taken into account. The non-ideal PSFs de- 
scribe the combined effect of the Broadband Filter Imager and 
the Solar Optical Telescope. In contrast to the ideal PSFs, they 
account for contributions caused by imperfections and contami- 
nation of the optical components in the telescope (e.g., impurity 
of lens material, dust on optical surfaces, microscopic scratches 
or micro-roughness) and stray light, which may originate, e.g., 
from reflection at baffles. The non-ideal PSFs are calculated by 
convolution of the ideal diffraction-limited PSFs and Voigt pro- 
files, which are chosen as model for the non-ideal contribution. 
For each wavelength channel, such theoretical point spread func- 
tions were produced for a grid of the parameters <x and y, which 
control the shape of the non-ideal Voigt-like contribution. These 
PSFs were then applied to artificial images of an eclipse and 
a Mercury transit. In a next step, the grid of resulting artificial 
intensity profiles across the terminators were compared to SOT 
observations of the Mercury transit from November 2006 and the 
solar eclipses from 2007. For each observed image in the blue, 
green, and red continuum channels of the BFI, the optimum PSF 
is indicated by the best-fitting artificial profile, resulting in the 
parameters [<x, y] and a usually small residual intensity offset. As 
the non-ideal contributions depend on a number of factors, the 
individual images produce a range of possible parameter combi- 
nations rather than a single solution. Nevertheless, the stray-light 
contributions are found to be best matched with Voigt functions 
with the parameters cr = 0'.'008 and y = 0'.'004, 0'.'005, and 
0'.'006 for the blue, green, and red continuum channels, respec- 
tively. 
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Fig. 2. Illustration of deconvolution using optimal (Wiener) fil- 
tering. The top panel shows a detail of a blue continuum disc 
center image. The left image is the flatfielded image, the middle 
image after deconvolution using the ideal PSF, and the right im- 
age after deconvolution using the non-ideal PSF from Paper I. 
All three images are scaled using the same intensity range: the 
minimum and maximum pixel values of the right panel. In the 
graph below, the thick solid line is the azimuthal average of the 
2D power spectrum of a blue continuum disc center image. The 
triple-dot-dashed line is the noise model; the dots represent the 
power spectrum with the noise model subtracted. The dash-dot 
line is the model for the signal in the absence of noise. Lines in 
medium thickness: the dashed line is the ideal MTF, the dotted 
line the non-ideal MTF. Both have their cut-off at the diffraction 
limit A/D = 0'.'19. Thin lines: the dashed line is the optimal fil- 
ter for the ideal PSF, the dotted line for the non-ideal PSF. Both 
lines relative to the y-axis scaling on the right side. 



2.3. Deconvolution 



The granulation images are deconvolved with (i) the ideal 
(diffraction-limited) PSFs and (ii) the non-ideal PSFs that are 
described in Sect. I2.2I We use all images at disc-centre for each 
wavelength channel (15 images each). The rectangular Hinode 
images were split into two square images for deconvolution in 
order to facilitate a straightforward execution and interpretation 
using standard FFT methods. After deconvolution, the average 
contrast of both sides is used for the analysis. 



The deconvolution was performed by means of an opti- 
mal ( Wiener) filter (see, e.g., Chapter 5.8 of lGonzalez & W oods 
2008). Figure[2]illustrates various aspects of the construction of 
the optimal filter. The graph shows the azimuthal average of the 
2D power spectal density (or "power spectrum", see Sect. [6]l of 
signal plus noise as the solid black line. In the signal-dominated 
regime, the power spectrum has a peak followed by decreasing 
power for the intermediate range of increasing wavenumbers k. 
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In the noise-dominated regime, the power increases with increas- 
ing wavenumber for spatial scales smaller than the diffraction 
limit AID < 0'.'19 (see the cut-off for the MTF). 

For the construction of the optimal filter, an assumption has 
to be made for the power spectrum of the sig nal in the absenc e 
of noise (see the discussion in Chapter 13.4 of Press et al. 1992). 
To this end, a noise model was constructed by extrapolating a 
fit to the noise tail into the signal region (triple-dot-dashed line). 
The subtraction of the noise model from the measured power 
spectrum is shown as black dots. The power spectrum displays 
an excess in the form of a less steep power fall-off around spatial 
scale 0'.'3. A similar excess is seen in all power spectra and be- 
comes more pronounced for the green and red continuum chan- 
nels. Since we do not observe such excess in power spectra from 
higher spatial resolution gran ulation images from the Swedish 1- 
m Solar Telescope (see, e.g.. lRouppe van der Voort et alj|2 004) 
and the Dunn Solar Telescope (see, e.g. JWoger e t al. 2008J), we 
decided to exclude this excess from the signal model. The ori- 
gin of the excess is unclear but we can exclude JPEG compres- 
sion artifacts (see Sect. 17.31 ). The signal model is derived from 
an extrapolation of the power spectrum to high wave numbers 
following the power law in the [0'.'6 - 1"] range (dash-dotted 
line). 

The azimuthal average of the optimal filter is also shown in 
the graph of Fig. [2] (thin lines, following the right y-axis scal- 
ing). The filter enhances (> 1) for spatial scales > 0'.'27 and 
suppresses noise for smaller scales. Although the shape of the 
optimal filter is significantly affected by in- or exclusion of the 
power excess in the signal model (including the excess shifts 
enhancement of the filter down to 0'.'2), the contrast is affected 
only little: 0. 1 % higher contrast] for the ideal PSF and 0.2% for 
the non-ideal PSF. This illustrates that the contribution of small 
spatial scales to the overall contrast is limited. That makes de- 
convolution by an optimal filter a robust procedure, where de- 
tails of the construction of the filter have only limited impact on 
the contrast values. 

While the deconvolution of disc-centre images produces 
rather robust results, the treatment of images at smaller /i is prob- 
lematic. In particular, images close to the solar limb cover a sig- 
nificant range in ju. Consequently, the properties of the stray light 
component vary significantly across the FOV. The same applies 
to the noise in power spectra that are used for the construction for 
the Wiener filter. It is therefore not possible to construct a Wiener 
filter that is optimal for the whole FOV. Retreating to deconvo- 
lution of smaller image sections with small intrinsic variation in 
ft would cause systematic errors due to neglect of the far PSF 
wings. The contrast values for deconvolved images away from 
the disc-centre would suffer a pronounced spread and possibly 
little significance only. We therefore restrict the deconvolution 
to disc-centre images. 

The three images in Fig. [2] show the effect of the decon- 
volution. The left panel shows the flatfielded image, the mid- 
dle panel after deconvolution considering the ideal PSF, and the 
right panel after considering the non-ideal PSF from paper I. All 
images are shown using the same intensity scaling. The effect of 
noise suppression for the two right panels is clear. The contrast 
is enhanced from 12.6 % for the flatfielded image to 17.9 % for 
the image deconvolved with the ideal PSF, and to 26.5 % after 
deconvolution with the non-ideal PSF. See Table [3] for statistics 
for all disc-centre images. 



3. Synthetic intensity 

3.1. Numerical simulations 

An overview over the models used here is given in Table [T] 
Most of our analysis will be based on three sn apshots of a re - 
cent 3D radiation hydrodynamic simulation by Steffen (2007), 
which was computed with CQ 5 BOLD ([ Frevtag e t al.l l2002t 
IWedemever et alj|2004t iFrevtag et alj|2008f) . The model which 
is identified as ChO hereafter, has a horizontal grid spacing of 
28 km and an extent of 11.2Mmx 11.2 Mm. It reaches from ~ 
-2400 km below the average height of optical depth Tsoo„ m = 1 
to +750 km above. The simulation code CO s BOLD employs a 
Riemann-type solver. Note that the effective spatial resolution 
of models produced with such a code is higher than for mod- 
els calculated with a finite-differences solver at same grid spac- 
ing. The radiative transfer in the simulation uses a frequency- 
dependent ("non-grey") long characteristics scheme with multi- 
group opacities with five bins. The latter are constru cted from 
updated MARCS opacities (see |Gustafsson et al. 2008, and ref- 
erences therein). 

For comparison, additional sets of synthetic images 
are calculated for CO 5 BOLD models and simulations by 
IStein & Nordlundl dl998l) . All sets consist of three snapshots 
each. The average contrast and standard variation for disc- 
centre images are given in Table [3] The abbr eviations Chi and 
Ch2 s tand for the non-magnetic models by Wedemever et al.l 
(2004) with frequency-independent ("grey") and frequency- 
dependent ("non-grey") radiative transfer, respectively. Cm is 
also a CQ 5 BOL D model but with a weak magnetic field of 
(B z ) = 10 G (ISchaffenberger et al.l [2006b . The models by 
Stein & Nordlund only differ such that SNh has no magnetic 
field, whereas SNm has (B~) = 250 G. A snapsho t of th e 



latter has already be en used by, e.g., Carlsson et a l. (2004) 



i Leenaarts et al. (2006), and bv lUitenbroek et al.l d2007l) . See also 
lstein& Nordlundl (12006). 



3.2. Intensity synthesis 

We use the spectrum synthesis code Linfor3D to 
produce synthetic intensity maps for the radiation 
(magneto-)hydrodynamic simulations listed in Table [T] 
Linfor3D (see http://www.aip.de/~mst/linfor3D_main.html) is 
originally based on the Kiel code LINFOR/LINLTE. It solves 
the detailed radiative transfer in full 3D for an input atmosphere 



Table 1. Numerical simulations used in this study. 



ID 



reference 



RF 



CCBOLD simulations 



parameters 
(B z ) b hor. ext.' 
[G] [Mm] 



ChO 
Chi 
Ch2 
Cm 



Steffen (2007) non-gr. 11.2x11.2 

Wedemeyer et al. (2004) grey 5.6 x 5.6 

Wedemeyer et al. (2004) non-gr. 5.6x5.6 

Schaffenberger et al. (2006) grey 10 4.8x4.8 



Stein & Nordlund simulations 



SNh 
SNm 



Stein & Nordlund (2006) non-gr. 6.3x6.3 
Stein & Nordlund (2006) non-gr. 250 6.3x6.3 



1 All contrast values and differences between them are given in per- 
cent points, notified by % (see Sect. 15.2b . 



a For each model it is indicated if frequency-independent ("grey") or 

-dependent ("non-grey") radiative transfer (RT) was considered. 
b magnetic field (B z ) 

c horizontal extent of the computational domain 
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Fig. 3. Comparison of synthetic and observed 
images for the blue continuum channel. Top 
left: observed deconvolved image, top right: 
synthetic image rescaled to the same pixelsize 
as the observed image; bottom left: same as 
top left but with contours enclosing the pix- 
els, which are brighter than in the synthetic im- 
age (//</> > 1.66). The image is made fainter 
for better visibility of the contours. Bottom 
right: histograms for the intensity distributions: 
original observed (dotted), deconvolved (dot- 
dashed), and original synthetic intensity (thick 
solid). 



under the assumption of local thermodynamic equilibrium 
(LTE). The computations are performed on a numerical grid 
that has a higher resolution compared to the input model. For 
consistency reasons, the same opacities and element abundances 
are used in the numerical simulations and in Linfor3D (e.g., 
updated MARCS opacities for model ChO, see Sect. 13. U . 

The green and in particular the red continuum band contain 
only a few mostly weak spectral lines so that their relative in- 
fluence on the radiation integrated over the whole wavelength 
band is negligible. The blue continuum channel includes a spec- 
tral line of neutral iron (Fe I) at 450.483 nm, which we explicitly 
take into account. The l ine parameters are taken from the Vienna 
Atomic Line Database dKupka et al.ll 1 9991) . The influence of the 
line on the radiation integrated over the full blue continuum filter 
turns out to be negligible. The change in contrast would be well 
below 0. 1 percent point with respect to pure continuum images 
and also the intensity histograms are barely affected. 

3.3. Degradation of synthetic maps 

The PSFs introduced in Sect. 12.21 are used to simulate the im- 
age degradation caused by the SOT/BFI optics. Each synthetic 
intensity map is convolved with the non-ideal PSF of the corre- 
sponding wavelength. For the model ChO, we also use the ideal 
PSFs and two different FOVs (see discussion in Sect. 17.31 ). 



4. Comparison of observed and synthetic images 

We begin our comparison of the Hinode filtergrams and the nu- 
merical simulations on the basis of disc-centre images. First, the 
de-convolved observations are opposed to the original synthetic 
intensity maps from model ChO in Sect. 14.11 The original ob- 



servations are then checked against the corresponding degraded 
synthetic maps in Sect. 14.21 

4. 1 . Synthetic images compared to deconvolved 
observations 

The synthetic images have a lower contrast than the PSF de- 
convoluted SOT/BFI images but are still within the 1-cr varia- 
tion of the observed values in the blue and green continuum. 
The differences are 1.7 %, 1.3 %, and 1.9 % for the three chan- 
nels, respectively. The comparison is illustrated in Fig.[3]for blue 
continuum images at disc-centre. 

In order to make a fair comparison, the synthetic image is 
filtered with a low-pass filter. The filtering accounts for the fact 
that the observations are taken with a telescope with finite aper- 
ture. Consequently, there can be no higher spatial frequencies 
present than given by the size of the aperture. In addition, the ob- 
servations contain noise at a level that affects the highest spatial 
frequencies (the power spectrum of the image is dominated by 
noise at frequencies higher than the diffraction limit). The low- 
pass filter was constructed in a way such that the power spectrum 
of the filtered synthetic image matches the power spectrum of the 
deconvolved image. After convolution, the synthetic image was 
re-sampled to the Hinode pixel scale. The observed image was 
selected to be as quiet as possible, selecting a region with as few 
magnetic bright points as possible. The top left image in Fig. [3] 
shows the deconvolved observed image, the top right image the 
filtered synthetic image. Both are shown using the same intensity 
scaling. In a visual comparison of the two images, one can dis- 
tinguish the observed from the synthetic image by the presence 
of a noise pattern that remains after application of the optimal 
filter. Other than this noise pattern, the two images are qualita- 
tively remarkably similar. 
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Fig. 4. Comparison of synthetic and observed 
images for the green continuum channel. Top 
left: original observed image, top right: de- 
graded synthetic image rescaled to the same 
pixelsize as the observed image; bottom left: 
original synthetic image; bottom right: his- 
tograms for the observed intensity (dotted), de- 
graded synthetic intensity (thin solid), and orig- 
inal synthetic intensity (thick solid). 



The lower right panel of Fig. [3] shows the histograms of 
the intensity values of the two images. In addition, the his- 
togram of the flatfielded observed image is shown (dashed line). 
Deconvolution by the non-ideal PSF has a pronounced effect on 
the intensity histogram. The histograms of the deconvolved ob- 
served and synthetic images are very similar. It should be noted 
that the low-pass filtering of the synthetic image has only limited 
effect on the histogram. Both the synthetic images and the de- 
convolved observations exhibit an essentially "double-peaked" 
distribution, although the peak at // (I) < 1 is much more pro- 
nounced than the secondary peak at //(/) > 1. The latter can 
easily be obscured and only be noticable as asymmetry in the 
distribution (see Fig . [6j> . The distribution can actually be decom- 
posed in a bright granular and a dark intergranular component 
for the simulations, when one considers the sign of the flow 
velocity (downflow/upflow). T his property has been kno wn for 
some time (see, e.g., Fig. 3 in IStein & Nor dlund 2000l). A re- 
markable difference between the two histograms, however, can 
be seen at high intensity values. These high pixel values (160 in 
total, or 0.2%) can be found as bright regions inside granules, 
enclosed by the contours in the lower left panel of Fig. [3] 



4.2. Degraded images compared to original observations 

Now we reverse the procedure by artificially degrading the syn- 
thetic images by convolution with the non-ideal PSFs described 
in Paper I. Using the best-fit PSFs with Voigt contributions, re- 
duces the contrast values of the synthetic image to 13.0 %, 8.6 %, 
and 6.0 % (average over three snapshots) for the blue, green, and 
red continuum, respectively (see Table [3). The resulting differ- 
ences between the average observational and degraded synthetic 
values are then 0.2, 0.3, and -0.2 percent points for blue, green, 



and red continuum, respectively. The above mentioned degrada- 
tion was performed with a PSF extending with a size given by 
the computational box of the model, i.e., 15 "A x 15 '.'4. The far 
wings of the PSF are this way neglected. We therefore repeat 
the degradation with a PSF with a size of 1 10" x 1 10", which is 
comparable to the FOV of the observations. The synthetic image 
is repeated periodically to fill the FOV. This procedure is war- 
ranted as periodic boundaries were used during the simulation 
run. The large synthetic images are now degraded with the full 
non-ideal PSF, including the far wings. At disc-centre, we find 
mean contrast values of 11.3 %, 7.4%, and 5.1 % for the blue, 
green, and red continuum, respectively (see Table [3]). This cor- 
responds to an effective reduction of the contrast ranging from 
0.87 to 0.85 of the value for the small FOV. On the other hand, 
uncertainties in the parameters for the model of the non-ideal 
PSF contributions have a tendency to increase the contrast again 
towards the value found for the small FOV (see Sect. 17.31 ). We 
conclude that the observed and synthetic snapshots agree well 
in terms of intensity contrast when considering the uncertainties 
inherent to the applied PSF. 

Figure|4]shows examples for an original observed and a syn- 
thetic image in the green continuum channel. The intensity dis- 
tributions of the original observed and the degraded synthetic 
image agree well. The distributions of individual images and 
subfields, however, vary to some degree in accordance to the 
variation in intensity contrast. Some cases show a narrower dis- 
tribution, some exhibit a number of particularly bright pixels (as 
in Fig.[3]for the deconvolved image), while the secondary com- 
ponent of brighter than average pixels (see see Sect. l4.U is some- 
times only noticeable as a distribution asymmetry. The distribu- 
tions are generally close to those found for the synthetic images. 
See Fig. [5] for a comparison of simulation and observation that 
includes all available disc-centre images. 
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Fig. 5. Comparison of the intensity distribution in all synthetic and observed disc-centre images for the blue (left), green (middle), 
and red continuum channel (right). Shown are the distributions for the original observed images (dotted), the deconvolved images 
(dot-dashed), the degraded synthetic images (thin solid), and the original synthetic intensity maps (thick solid). 



5. Centre-to-limb variation of the continuum 
intensity distribution 

In this section, the properties of the observed and simulated in- 
tensity distribution and contrast are investigated as function of 
heliocentric position. 

5. 7 . Observed intensity histograms. 

The intensity distribution of the original (reduced) SOT images 
is illustrated in Fig. [6] for the three continuum wavelengths. The 
pixels of all images are sorted by disc position fi into bins of 
size Afi = 0.05. The upper row of the figure shows the his- 
tograms for all bins, whereas the rows below show histograms 
for selected heliocentric positions. The intensity distribution be- 
comes narrower with increasing wavelength and with decreasing 
fi. The histograms are relatively symmetric close to the limb but 
exhibit a growing asymmetry towards disc-centre. The asymme- 
try is most pronounced at [i = 1 .0. Generally, there seems to be 
an increasing number of pixels with intensities brighter than av- 
erage, which cause the distribution to be broader at disc-centre. 
Consequently, the peaks of the histograms move from around 
1/ < I >~ 1.0 at the limb to smaller values at disc-centre. 

Removing the PSF (see Sect. I2.31 l from the images has a 
dramatic effect on the resulting intensity distribution. It broad- 
ens, covering a much larger range in //(/), and exhibits a subtle 
secondary peak for brighter than average pixels. More details 
are discussed in comparison to synthetic images in Sect. [4] and 
Fig. [3] 

5.2. Observed intensity contrast 

Definition. The intensity contrast is defined as 

V*Z J|y (/(*,3» f 

O/rms = — , (1) 

V )x,y 

where I(x,y) are the intensity values of the individual pixels 
with the average {I) x , y over all N pixels. 

Contrast of reduced images. The contrast values are listed for 
the three considered BFI continuum bands and different helio- 
centric positions fi in Table [2] The average contrast values for 
the SOT filtergrams are derived from regions of equal fi with a 



bin size of Afi = 0.05 for all images. At disc-centre, the im- 
ages are within the same //-bin, whereas image segments are 
sorted into the bins closer to the limb. The complete images 
at disc-centre (15 for each wavelength), produce the follow- 
ing average contrasts and standard deviations: (12.8 ± 0.5)%, 
(8.3 + 0.4) %, and (6.2 ± 0.2) % for blue, green, and red con- 
tinuum, respectively. The intensity contrast strongly decreases 
towards the longer wavelengths. The same behavior is found 
for all fi. For each wavelength, the contrast also decreases from 
centre to limb, i.e. with decreasing fi, which can also be seen 
from the width of the histograms discussed in Sect. 15.11 This 
centre-to-limb variation of the contrast is shown in Fig. [7] which 
can be interpreted as grey-shaded histograms. The histograms 
for the individual bins exhibit a maximum close to the corre- 
sponding average contrast (solid circles, cf. Tabled. The distri- 
bution, however, cannot be described by means of a, e.g., simple 
Gaussian. The contrast distribution for the blue continuum even 
shows a small number of cases with enhanced contrast, which 
manifest as a secondary ridge in Fig.|7^. This is most likely con- 
nected to the influence of magnetic fields, e.g., in the form of 
network bright points, which are unavoidable for the large size 
of the used FOV, even though enhanced network regions were 
explicitly excluded from the sample (see Sect. [4]). Small errors in 
the determination of the limb and also the limited spatial resolu- 
tion of SOT introduce uncertainties for the determination of the 
contrast very close to the limb (see Sect. 17.3b . Data for fi < 0.07 
should thus be interpreted with caution. 

Contrast after PSF deconvolution. Removing the instrumen- 
tal influence by deconvolution with the PSF has a dramatic ef- 
fect on the contrast. The results are summarised in Table [3] The 
contrast of blue continuum images at disc-centre is increased 
by 5.2 percent points after deconvolution with the diffraction- 
limited ideal PSF and another 8.7 points, when the full non-ideal 
PSF is used. In total, the original contrast is more than dou- 
bled to (26.7 + 1.3) %. The green and red channel show sim- 
ilarly large effects. The contrast increases from (8.3 + 0.4)% 
to (19.4 + 1.4)% in the green and from (6.2 + 0.2)% to 
16.6 + 0.7) % in the red channel. 

Error estimate. It is not trivial to give a reliable error estimate 
for the empirical contrast values. Next to systematic errors due 
to instrumental properties (PSF, straylight, see Sect. I7.3K there 
are statistical errors arising from the selection of the pixel en- 
sembles. The contrast values given in Table [2] are derived from 
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Fig. 7. Centre-to-limb variation of the intensity contrast for the blue (left column), green (middle), and red continuum (right). The 
contrast values for the SOT wide-band filtergrams are derived from regions of equal fi (bin size Afi = 0.05) for all original images. 
The individual histograms are combined to a set of /i-dependent histograms ("density map"), which are shown as grey scales. The 
maxima of the histograms are marked with solid circles for each fi bin. For comparison, the contrast values for the 3 ChO snapshots 
are shown as double triangles, while the solid line represents the mean centre-to-limb variation of model ChO. 



pixels with similar heliocentric position fi, effectively forming 
segments of circles in the images. At disc-centre, the change in 
fi over one BFI image is so small that the complete images in- 
stead of segments are used. The rms values in Table [3] thus rep- 
resent the variation among whole images. The smaller the fi bin- 
size, the smaller the number of the included pixels over which 
intensity mean and variation are computed. This can lead to a 
small increase of the scatter in contrast at small fi, while the re- 
sults for positions closer to the disc-centre are less susceptible. 
Another way to determine a statistical error is to compare the 
variation of the contrast in square sub fields. This error increases 
slightly with decreasing subfield size. While the same rms vari- 
ations than in Table [3] are found for large field sizes, the spread 
increases slightly to 0.7 %, 0.6 %, and 0.3 % for blue, green, and 
red continuum, resp., when dividing the images into fields of 
only 10 "x 10". 

For the construction of the optimal filter used for deconvo- 
lution, assumptions have to be made for the signal model. As 
already stated in Sect. 12.31 the two different signal models con- 
sidered here result in contrasts that differ by only ~ 0.2 percent 
points. The deconvolution procedure seems to be rather robust. 
We degrade the deconvolved images again with the same PSF in 
exactly the way as it is done for the simulations in Sect. [4] The 
original observed images and the back-transformed images have 
contrast values that differ by much less than 0.1 percent points. 



5.3. Synthetic intensity contrast 

For fi — 1, the three selected snapshots of the model ChO give 
an average contrast of 25.0 %, 18.1 %, and 13.8 %, at the central 
wavelengths of the blue, green, and red continuum channels of 
the BFI, respectively. The variation among the individual snap- 
shots is for all wavelengths of the order of 0. 1 percent point. 

Generally, there is a spread of a few tenths of a percent point 
among individual snapshots. All non-magnetic CO s BOLD mod- 
els and also SNh produce very similar contrasts. Obviously, the 
model ChO analysed in more detail in Sect.|4]is no exception but 
representative for state-of-the-art granulation simulations. Also 
the weak field model Cm shows no systematic difference with 



respect to the non-magnetic models. In contrast, the model with 
the high field strength (SNm) exhibits systematically lower con- 
trasts with respect to SNh at disc-centre. The values nevertheless 
remain roughly of the same order as for the other models but be- 
come significantly larger at smaller fi (see Fig. [8]). 

The contrast of all non-magnetic and weak field models de- 
creases monotonically with decreasing fi. The difference be- 
tween the individual simulations can essentially be described as 
a small contrast offset. An interesting exception is the higher 
field strength model SNm. Starting with a lower value than SNh 
at disc-centre, the contrast grows with decreasing fi until a peak 
is reached at fi — 0.3. For fi < 0.8, SNm exceeds the contrast 
of SNh. The blue continuum band of SOT in Fig. |7h indeed 
shows a secondary group of image segments with a contrast ex- 
cess, which starts to divert from the main distribution already 
at disc-centre. It is not obvious at which fi this subset of SOT 
observations reaches maximum contrast but it could well be for 
smaller ft than seen in the model SNm. The latter represents a 
plage region with (\B\) = 250 G, whereas the effect is not clearly 
discernible for the weak-field model Cm. We therefore argue 
that the subset of SOT data with contrast excess might contain 
a small number of magnetic field structures with smaller field 
strengths. 

The deviation of the contrast is connected to the appear- 
ance of solar faculae close to the limb. Faculae can be under- 
stood as a line-of-sight effect caused by magnetic elements in 
the intergranular lan es. This topic has be e n discussed extensivel y 
in recent stud i es bvlKeller et all d2004l) . ICarlsson et"ai1 ( 12004), 
ISteinerl (l2005l).lDe Pontieu et alJd2006l) . andlBerger et alJd2 007 ; 

1 C0051) . ISanchez Cuberes et al 



but als o by Iffirz berger & WiehrJ (120051) . ISanchez Cube 
(|2003[) . ISutterlin et al l (1 1999b . and many more. See |0: 



rtiz et al 



(2002) for an empirical investigation of the centre-to-limb vari- 
ation of the intensity contrast for solar photospheric faculae. 

The additional high-contrast ridge seen in the blue is not vis- 
ible for the other two SOT channels in Fig. [7] This might partly 
be a selection effect as we explicitly selected quiet Sun images, 
avoiding enhanced magnetic field whenever possible. The im- 
ages were taken during solar minimum so that the Sun was very 
quiet and not many pronounced active regions were present. Still 
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Fig. 8. Centre-to-limb variation of the intensity contrast for all models: top: original synthetic images, bottom: after degradation with 
the non-ideal PSFs. The columns show the results for the blue (left), green (middle), and red continuum band (right), respectively. 
The contrast values for the original SOT wide-band filtergrams are represented as gray-shades for reference (cf. Fig. [7j in the lower 
row. The snapshots of each model are represented by symbols, while the lines show the mean centre-to-limb variation. Please note 
that only the actual simulation box sizes are considered for the degradation of the synthetic images shown here. Expanding it to a 
FOV comparable to the SOT observations would lower the individual contrast values but would give practical problems for low yU 
(see Sect. [73), 



one has to consider weak magnetic fields on small spatial scales 
in the quiet Sun. 



The small contrast differences among the non-magnetic 
models are partly systematic effects due to differences 
in the radiative transfer schemes a nd in the opacity ta- 
bles, e.g., PHOEN IX/OPAL opacities (lHauschildt etal.lll997t 
Igles ias et al.l lT992) for Chi . On the other hand, already the ran- 
dom selection of three snapshots influences the mean contrast 
by the order of a few tenths of a percent point. A larger set of 
63 images for Ch2, for instance, produces contrasts in the blue 
continuum between 23.7 % and 26.7 % with a mean of 25.2 % 
and an rms variation of 0.7 %. While the mean is reduced by 0.3 
percent points for the larger sample, the rms variation increases 
from 0.5 % to 0.7 % (cf. Table [3). A determination of the mean 
granulation contrast with a precision down to the order of a few 
0.1 % obviously requires a statistically significant sample of im- 
ages rather than a single snapshot. 



6. Spatial power spectral density 

6. 1 . Observed power spectral density 

The power spectral density Pk (PSD or short "power spectrum") 
is calculated by integrating the amplitude square H 2 of the two- 
dimensional Fourier transform of an image over rings of constant 
wavenumber k. We use wavenumbers defined as 



k= yjk 2 . + k 2 -2n tJl x ~ 2 + L y - 



(2) 



where L x and L y are the spatial wavelengths in x- and y-direction. 
Please note that here the one-sided PSD is used, i.e. the Fourier 
amplitudes H for both the ne gative and the po sitive wavenum- 
bers are considered (see, e.g. JPress et al.|[l992h . As the Parseval 
theorem is fulfilled, the total power for the intensity // (/) is re- 
lated to the contrast of the original image: 



P k dk' = (£/ rms ) 2 



(3) 



Power spectra are calculated for granulation images at disc- 
centre for all three considered SOT continuum bands. The az- 
imuthally integrated power spectra are plotted versus k in the 
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Fig. 9. Power spectral density P (left column) 
and kP k (right column) for //(/) as functions of 
spatial wavenumber k for the three BFI contin- 
uum wavelengths at disc-centre: blue (top row), 
green (middle), and red continuum (bottom). 
The observational data (thin solid) is compared 
to the the degraded synthetic data (thin dot- 
dashed), while the original synthetic images 
(thick dot-dashed) are checked against the ob- 
servations after deconvolution with the ideal 
PSFs (thick dotted) and the non-ideal PSFs 
(thick solid). The slopes of the power spec- 
tral density can be approximated with power- 
laws of the form oc kf p . The exponents are 
Sp ~ 1.0 for small k and between ~ -4.0 and 
~ -5.0 for large k (dotted lines). For kP k , the 
exponents are correspondingly s kP = s P + 1 
(see text for details). The dotted lines repre- 
sent the power law fits to the observational data. 
The solid vertical lines indicate the wavenum- 
ber corresponding to the diffraction limit of the 
telescope. 



left column of Fig. [9] The curves show power maxima at a 
spatial scale of ~ 1.7 Mm, 1.8 Mm, and 1.9 Mm for the blue, 
green, and red continuum, respectively. That scale is somewhat 
larger than the most frequent granule diameter, which is typi- 
cally of the order ~ 1.4 Mm in the SOT images analysed here. 
The latter is in very good agreem ent with earlier resu l ts suc h 
as, e.g., with the 1 '.'9 reported by iBrav & Loughhead] dl977l) . 
The spectra also show signs of enhanced power density at spa- 
tial scales > 20 Mm that cor r espo nd to super-granulation. This 
is in line with Rieu tord et al.l d2008l) . who derive the kinetic en- 
ergy density spectrum by granule tracking. Except for these ad- 
ditional contributions, the slopes can roughly be approximated 
with power laws of the form P k oc k Sp . The exponents sp for 
small k < 2.0 Mm 1 are of the order of 1.0, i.e. the power in- 
creases essentially linearly with wavenumber. Fits of the ob- 
servational data produce exponents that decrease slightly with 
wavelength from 1.0 for the blue, to 0.9 for the green, and 0.8 
for the red channel. Power law fits are also possible for large 
wavenumbers although the exponents change with k and thus 
depend on the considered wavenumber range. In the range from 
9Mirr' to 18Mtrr', we find exponents of roughly -5.0, -4.5, 
and -3.7 for blue, green, and red continuum, respectively. The 
slopes get steeper at even higher wavenumbers before the power 
spectral density is dominated by noise and increases again. The 
contributions to the power spectral density at these wavenumbers 
are, however, of minor importance only. 

The right column of Fig. [9] shows the quantity kPk- 
According to Eq. (0), the area below the curve is equal to the 
intensity contrast squared. The distributions have clear peaks 
corresponding to those in the power spectra in the left column 
but at somewhat smaller spatial scales: 1.4 Mm for blue and red 
continuum and at 1 .5 Mm for the green channel. These scales co- 



incide with the most frequent granule diameter of ~ 1 .4 Mm (see 
above). The slopes can again be approximated with power laws 
of the form kP k oc k Spk . The exponents sp k should correspond to 
s kP = sp + 1 and are indeed found to be in the range between 
-4.0 and -2.7. 

6.2. Synthetic power spectral density 

Not only the intensity distribution but also the power spectral 
density of model ChO is in line with other simulation s as can 
be seen, e.g., from Fig. 3 bv lStein & Nordlundl J2000h . It illus- 
trates that model ChO is representative of state-of-art numerical 
simulations of solar granulation. 

6.3. Comparison of observed and synthetic PSD 

Figure [9] demonstrates that the power spectral density of the de- 
graded synthetic images matches well the observational data de- 
scribed in Sect. |6]for all three continuum bands. The power law 
fits for the degraded synthetic data produce exponents that agree 
very closely with the empirical exponents. The slope for larger 
wavenumbers in the red continuum, however, is better repre- 
sented with an exponent of sp w -4.0, which is slightly smaller 
than the observational analogue. At these wavenumbers, the in- 
strumental image degradation becomes significant and, by sup- 
pressing power at these scales, influences the slopes. The result- 
ing steep decrease in power spectral density is also present for 
the degraded synthetic images. 

The power spectral density for the original synthetic images 
is generally higher corresponding to a higher intensity contrast 
(see Eq. (01). The power law exponents are determined as 1.0 
for small wavenumber and 4.5 at large wavenumbers for all 
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Fig. 6. Intensity histograms for the blue (left column), green 
(middle column), and red continuum (right column) as function 
of disc-position p. Top row: Histograms for all p with the rel- 
ative occurrence as color/grey-shades. Rows below: Histograms 
for p. = 0.1, 0.3, 0.5, 0.7, and 1.0, respectively. The lines mark 
//(/) = 1 (dot-dashed), the histogram peaks (dashed), and - in 
the lower rows - the histograms mirrored at the histogram peaks 
(dotted). The latter reveals the asymmetry of the distribution. 



three wavelengths and thus still agree with the observations. The 
slight wavelength dependence found for the observed and the de- 
graded data is not present or at least less pronounced in the orig- 
inal synthetic intensity maps. Furthermore, the power spectral 
density of the latter continues to decrease further with wavenum- 
ber well into the domain where the observational signal (and also 
the degraded synthetic one) is drowned in noise (see Fig. [9j, i.e., 
k > 40 Mm" . Due to that steep slope, the integral of the PSD is 
changed only slightly when going towards higher wavenumbers. 
According to the normalisation in Eq. (|3), a small change of the 
integral corresponds to a small change in intensity contrast. The 
contribution to the intensity contrast is therefore very small at 
high wavenumbers but large for granulation scales correspond- 
ing to smaller k. 

We illustrate this with an experiment in which we artifi- 
cially filter away the small scales in synthetic images for the blue 
band. It is done by Fourier-transforming an image, eliminating 
all Fourier-components for wavenumbers larger than an imposed 
cutoff £ cuto ff, and back-transforming. The power spectral density 
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Fig. 10. Fourier-filtering of synthetic images at 450 nm. Intensity 
contrast 6I Tmii as function of the imposed spatial cutoff cr cutoff (a) 
and its derivative (b). The power spectral density (c) and the in- 
tensity distributions (d) are plotted for the original images (dot- 
ted lines) and two selected cutoffs of 0'.'3 (dot-dashed) and 076 
(solid), which are marked as open circles and/or vertical dotted 
lines in panels a-c. 



of the resulting image is thus truncated at k culo ff (see Fig. [TUb). 
The intensity contrast of the filtered images is plotted in Fig.fTOk 
as function of the spatial cutoff scale cr cuto ff = 27r/A: cu toff ■ The 
contrast decreases monotonically with the cutoff. Towards the 
smallest scales, 5I rms saturates to the value of unfiltered images 
(25.0%). The asymptotic approach to the final value already 
starts at scales that are well resolved in the numerical model (grid 



Table 2. Intensity contrast (5/ rms in percent points for the original 
observations in the three considered BFI continuum bands. 
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(7.5) 


5.5 


(5.6) 


0.80 - 0.85 


11.9 


( 11.6) 


7.7 


(7.4) 


5.7 


(5.6) 


0.85 - 0.90 


12.1 


( 12.2) 


7.9 


(7.8) 


5.9 


(5.8) 


0.90 - 0.95 


12.4 


( 12.5) 


7.9 


(7.8) 


6.1 


(6.2) 


0.95- 1.00 


12.8 


( 13.1) 


8.3 


(8.4) 


6.2 


(6.2) 



a The average contrast in each bin is given next to the peak of the cor- 
responding contrast histograms (in parentheses, compare to Fig.[7j- 

b heliocentric positions p; The p bins (Ap = 0.05) are the same as in 
Fig.H] 
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size ~ 0'.'039). The change of the slope can be seen from deriva- 
tive 5/rms/crcutoff in Fig. [TUb . (The derivative at larger cr cuto ff is 
noisy due to the numerical treatment of the discrete Fourier fil- 
tering.) The biggest change of the contrast occurs at scales be- 
tween > 0'.'5 and ~ 2 '.'2, which coincides with the "granulation 
peak" in the PSD. In contrast, the change towards the small- 
est scales decreases quickly. Removing the smallest scales has 
apparently only a small effect in terms of contrast. A cutoff of 
o"cutoff = 0"5 still produces a contrast of 24.2 %, i.e., a reduction 
of only 0.8 percent points. The absence of small scales is more 
obvious when looking at the intensity histograms in Fig. [TOt l. 
The distributions for the unfiltered images and for <x cutoi f = 0'.'3 
are very similar, whereas cr cuto ff = 0'.'6 already produces no- 
ticeable differences. The otherwise clearly visible double-peak 
transforms into a single-peaked asymmetric distribution. The re- 
duced occurrence of the brightest and darkest pixels implies that 
they are connected to the smallest scales. The darkest pixels, 
which constitute the 'major component of the intensity distribu- 
tion (at / < (/), cf. Sect. 14. H and Figs. [6] and |3j, are found in the 
narrow intergranular lanes. In the case of SOT, the lanes seem 
to be sufficiently resolved. We therefore expect no change of the 
contrast numbers for deconvolved ima ges from the use of tele- 
scopes with larger apertures. See, e.g.. iRao et al.l d200ll) for an 
analysis of the granulation contrast as function of the telescope 
aperture. 

Finally, we compare the PSD derived from the deconvolved 
observations with those for the original synthetic images, which 
are all plotted as thick lines in Fig. [9] The dotted and solid lines 
represent the PSD for deconvolution with the ideal and the non- 
ideal PSFs, respectively. The original simulations are marked 
with the dot-dashed lines. While deconvolution with the ideal 
PSF already produces a significant increase in power througout 
the observationally accessible wavenumber range, the non-ideal 



Table 3. Average contrast <5/ lms and standard deviation of the 
original and the deconvolved disc-centre observations in com- 
parison with the synthetic images. 



Sin 



[%] 


450. 


45 nm 


555.00 nm 


668.40 nm 


observation 


original 


12.8 


±0.5 


8.3 ± 


0.4 


6.2 


±0.2 


ideal PSF 


18.0 


±0.8 


12.3 ± 


0.8 


9.9 


±0.4 


full PSF 


26.7 


± 1.3 


19.4 ± 


1.4 


16.6 


±0.7 


simulation ChO (Steffen 2007) a 


original 


25.0 


± 0.1 


18.1 ± 


0.1 


13.8 


±0.1 


ideal PSF 


16.7 


±0.3 


11.4 ± 


0.2 


8.2 


±0.2 


non-ideal PSF (1) 


13.0 


±0.3 


8.6 ± 


0.2 


6.0 


±0.2 


non-ideal PSF (2) 


11.3 


±0.2 


7.4 ± 


0.2 


5.1 


±0.2 




other simulations (original) 


J 






Chi 


25.8 


±0.5 


19.0 ± 


0.5 


14.6 


±0.5 


Ch2 


25.5 


±0.5 


18.4 ± 


0.5 


13.9 


±0.5 


Cm 


25.4 


± 0.1 


19.0 ± 


0.5 


14.7 


±0.1 


SNh 


25.2 


±0.6 


18.4 ± 


0.4 


14.0 


±0.3 


SNm 


22.4 


±0.5 


16.0 ± 


0.3 


12.3 


±0.2 



The effect of the deconvolution with an ideal and a detailed PSF 
is specified for the model by Steffen (ChO). The extent of the de- 
tailed PSF corresponds to (1) only the size of the model and (2) a 
periodic repetition of the model to a size of 1 10" x 1 10", which is 
comparable to the FOV of the observations. 

For comparison, the original, i.e., non-degraded contrast values for 
sets of CO s BOLD models (Chi, Ch2, and Cm) and simulations by 
Stein & Nordlund (SNh and SNm) are given. 



PSFs bring the observations very close to the simulation results. 
The PSDs agree well for the blue and green continuum. For the 
red continuum, however, we remain with a small discrepancy, 
which is consistent with the differences in terms of intensity 
contrast (see Table. [3). The sharp decrease of the deconvolved 
PSDs close to the theoretical cutoff wavenumber (vertical lines 
in Fig. [Tol l results in a significant power deficit with respect to 
the simulations. This is simply due to the fact that there is no 
further information content in the images at the corresponding 
spatial scales due to the limited spatial resolution. Consequently, 
the power contributions at wavenumber close to the spatial cutoff 
are lost. 



7. Discussion 

7.1. Summary. 

Both simulations and observations are found to have a consid- 
erable variation in contrast. It is mostly a statistical variation 
between different snapshots / images. And still there are system- 
atic, physica l sources to be considered. lTitle et al.ldl989L see also 
ISpruit et al.l 1990) have shown that removing the intensity mod- 
ulation due to f- and p-mode oscillations in the data taken with 
SOUP at the former Swedish Vacuum Solar Telescope reduces 
the rms intensity variations by 35%. The relative contribution 
of oscillations to the contrast depends on the spatial resolution 
as oscillations are correlated with certain sp atial scales of whic h 
some possibly cannot be resolved (see, e.g. JSpruit et al.lll990b . 

We conclude that in view of the mutual variation ranges one 
cannot expect a much better match of synthetic and observa- 
tional continuum intensity contrast than has been achieved here. 

7.2. Comparison with other studies 

Observed contrast values for quiet Sun granulation vary signifi- 
cantly throughout the literature, indicating the large uncertainties 
arising from often unknown properties of the employed in stru- 
ment and seeing effects. See Sanc hez Cuberes et a il (120001) for 
an extensive survey of observed contrast values. Here we discuss 
some of the more recent studies that report granulation contrast 
me asurements. 

iKhomenko et alj d2003h determined a granulation contrast 
between 5 and 6 % for the continuum around 500 nm from ob- 
servations with the VTT Our SOT data rather indicate a contrast 
valu e of ~ 1 1 % for flat-fiel ded images at this wavelength. Note 
thatP^omenk o~et~aT] d2003l) make no corrections for instrumen- 
tal and seeing effects. 

Also us ing the VTT, de B oer (1995) found a contrast of 
14.5% (cf . iPerez Rodriguez & Kneerl |2002|) in speckle recon- 
structed images at 550 nm, which is significantly lower than the 
19.4 % we derived for the green continuum channel of SOT/BFI 
after deconvolution with the non-ideal PSF. 

Van Noort et al. d2005l) report granulation contrast numbers 
for SST observations in the continuum at 436.4 nm. Their im- 
ages are restored with the MOMFBD image restoration tech- 
nique and have maximum contrast values slightly below 11%. 
These values are significantly lower than the 12.8% that we 
find for the SOT/BFI blue continuum at slightly longer wave- 
length before deconvolution (see Table [3]). The PSFs they de- 
rived from MOMFBD processing are based on sub-images of 
size ~5"x5"and can therefore not account for the far wings 
of the PSF. No correction for scattered light is applied through 
MOMFBD processing. As we will discuss later in this section, 
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this is probably the main reason why their contrast values are 
rel atively low. 

lUitenbroek etail (120071) compare speckle-reconstructed G- 
band images of the quiet Sun taken with the Dunn Solar 
Telescope (DST) with a synthetic image calculated from a 
snapshot of radia t ion m agnetohydrodynamical simulation by 
IStein & Nordlundl (120061) with an initially horizontal magnetic 
field with <|B()|) = 30 G. They find a considerable discrep- 
ancy between the observed contrast of 14.1% for their best 
reconstructed i mages, and the synthetic contrast of 21.5%. 
Uiten broek et al. also analyse a synthetic image from the same 
Stein & Nordlundl simulation with (B z ) = 250 G that we use here 
(see Sect. [XT). They determine that the contrast of this snapshot, 
which corresponds to a plage region, is only 16.3% in the G- 
band. The fact that the effect of increasing the magnetic field 
leads to a decrease in the contrast, leads them to speculate that 
the quiet Sun might contain a larger amount of weak magnetic 
field than the 30 G assumed for the quiet Sun model. 

Our comparison of synthetic contrast values in the contin- 
uum around A = 450 nm in Sect. 15.31 also shows that a higher 
average magnetic field strength lowers the c ontrast (see Tab led 
but to much smaller extent than found by Uitenbroek et al.l in 
the G-band. This suggests that the presence of a large amount of 
spectral lines in the G-band is an important factor for the value 
of the granulation contrast as one compares magnetically active 
and quiet regions. 

We analyse speckle-recon structed images ob tained with the 
Dutch Open Telescope (DOT, Rutten et al. 2004), which we ob- 
tained from the DOT online data base. We find contrast values 
ranging between 14.4% to 18.6% for the G-band and values 
in the range from 14.9% to 19.8% for the DOT blue contin- 
uum at 432 nm. The code used to reconstruct the DOT images 
is bas ed on an IDL code first mentioned in Ide Boer & Kneerl 
(QUI. While some of the DOT contrast values are close to 
the DS T result from | Uitenbroek et al.l (120071) . others are clearly 
higher. lUitenbroek et alj find values around 14 % for G-Band im- 
ages obtain ed with the DST using the sp eckle code described 
in detail in IWoger & von der Liihel J2008). This code is differ- 
ent from the cod e used for the aforemen tioned DOT observa- 
tions. The code by Woge r & von der Lti he takes into account the 
usage of an Adaptice Optics (AO) system. The functions used 
for the Fourier amplitude c alibration are modified for Ada ptive 
Optics usage accor ding to I Woger & von derLuhd (120071) . bee 
Wog er etall (12008) for an analysis of the code's photometric ac- 
curacy in comparison to the Hinode satellite. 

The previous paragraphs illustrate that there is a considerable 
spread among the contrast values of reconstructed images due to 
inherent sources of uncertainty. In the case of speckle interfer- 
ometry, uncertainties in the Fourier amplitude calibration due to 
an incorrect estimation of the Fried parameter and/or (if appli- 
cable) the AO correction performance - and thus the application 
of the wrong calibration function - can lead to a large scatter in 
contrasts for the reconstructed images. 

In general, instrumental straylight is not accounted for in 
neither the speckle nor the MOMFBD reconstruction process, 
and needs to be treated separately. Therefore, the contrast values 
from post-facto reconstructed images should at best be compared 
to those disc-centre images of SOT that have been deconvolved 
with the ideal PSF only (e.g., (18.0 + 0.8) % for the blue contin- 
uum). But even then it is not guaranteed that the effect of seeing 
is completely removed while this source of uncertainty is of no 
concern for the space-bourne measurements with Hinode. 

The apparent contrast dilemma can be solved by reviewing 
the problematic treatment of non-ideal PSF contributions. The 



determination of the straylight level is very difficult, in particu- 
lar for ground-based observations. Consequently, details of the 
straylight properties are often unknown. A frequent strategy to 
account for straylight is to simply assume a constant contribution 
a, which is equally spread over the FOV. It corresponds to an in- 
tensity offset, which decreases the intensity contrast to (1 - a) 
of its original value. lUitenbroek et al.l (120071) conclude that they 
would need a straylight level of a = 0.34 in order to reproduce 
the contrast in their synthetic images. They consider this a factor 
3 too high for the DST. We also determine a from the difference 
between SOT images deconvolved with just the ideal diffraction- 
limited PSF and the full non-ideal PSF (see Table [3). In case 
of the blue continuum at disc-centre, the contrast difference of 
8.5 percent points translates into a — 0.33. For the green and red 
continuum, we find 0.37 and 0.41 , respectively. Alt hough these 
values are very close to the one by lUitenbroek et aU they should 
not be interpreted as actual straylight levels. The eclipse and 
Mercury transit observations, on which the SOT PSF with non- 
ideal Voigt-function contributions is based on, show residual in- 
tensities of the order of a few percent points only (see Paper I). 
Obviously, the assumption of a constant straylight contribution 
a is misleading. A constant offset would not only influence the 
PSF wings as intended but also the PSF core. For SOT, the stray- 
light level was actually shown to be dependent on the mean count 
rate, i.e., the overall light level in the instrument, and with that 
on the heliocentric position of the FOV. We therefore argue that 
a more detailed straylight model is crucial for detailed compari- 
ons of the granulation contrast and the intensity distributions in 
ge neral. 

iDanilovic et al.l (2008) obtained a continuum intensity map 
at 630 nm wavelength with the spectro-polarimeter (SP) of SOT. 
For this single map, they report a contrast of 7.0%, a value which 
is in line with our results from SOT/BFI. After detailed mod- 
elling of the SOT/SP optical syste m, they degrade a synthetic 
image from a MHD simulation by Vogler & Schiissler (2007) 
taking into account a detailed PSF. For the degraded synt hetic 
image, they find a contrast of 8.5%. IDanilovic et al.l (120081) ar- 
gue that a slight defocus of the SP instrument can account for a 
further decrease of the synthetic contrast of 1%. The remaining 
0.5% discrepancy is then attributed to straylight and imperfec- 
tions in the optical system, although one could also argue for 
an intrinsic spread in the observed contrast values. Our results 
point at a range of observed contrast values on the order of a few 
tenths of a percent around 630 nm. 

In their recent paper, Truiillo Bueno & Shchukina (2009) 
present the rms granulation contrast derived from detailed 
3D radiative transfe r for the 3D hydrodynamic model by 
Aspl und et al.l (120001) . Judging from their Fig. 9, the synthetic 
intensity maps have contrasts of roughly 24 %, 17 %, and 12 % 
for wavelengths corresponding to the blue, green and red chan- 
nels of the BFI. These values are systematically lower than for 
the synthetic maps considered here (see Table |3). The differ- 
ence, however, is well within the variation range of individual 
observed images and uncertainties discussed in Sect. 17.31 

7.3. Additional sources of uncertainty 

Spectrum synthesis: In Linfor3D, pure local thermodynamic 
equilibrium is assumed, i.e. scattering is not taken in to ac- 
count. In their Fig. 8, iTruiillo Bueno & Shchukina! (|2009) 
present the centre-to-limb variation of the continuum intensity 
at A = 300 nm. They demonstrate that accounting for the scatter- 
ing contribution to the source function removes the discrepancy 
between the observed and synthetic intensity. The discrepancy, 
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however, is only of the order of a few percent points for small 
/j and vanishes towards disc-centre. Scattering is even less im- 
portant for the wavelengths considered here. We therefore argue 
that the neglect of scattering in Linfor3D has only small if not 
negligible influence on our results. 

Point Spread Function. As the conditions may vary some- 
what between different observations, the corresponding PSF also 
varies. Consequently, each image would require its very own 
specific PSF for deconvolution. In practice, it was only possi- 
ble to determine PSFs for BFI/SOT that represent the best-fit 
cases for a large range in observational conditions (see Paper I). 
They may be considered as representative of a typical observa- 
tional situation. The analysis of the deconvolved images should 
therefore be based on a statistically significant set of images. 

Also the degradation of the synthetic images with the PSF 
bears sources of uncertainty. The spatial dimensions of the sim- 
ulations used here are smaller than for the observed images and 
the extent of the PSFs from Paper I. The convolution of a syn- 
thetic map with a PSF therefore mainly takes into account the 
central parts of the PSF while the far wings are effectively ne- 
glected. As mentioned already in Sect. 14.21 one should consider a 
FOV that is comparable to the one of the observation. Periodical 
extension of the synthetic maps to a size of 110" x 110" and 
degradation by a correspondingly larger PSF, lead to a reduced 
intensity contrast for the model by Steffen (see TableO. A lower 
contrast of the synthetic maps seems actually to be in line with 
the difference between the contrasts of the deconvolved observa- 
tions and the original synthetic maps. 

The influence of a too small FOV can also be seen for 
the other models in Fig. [8] Except for the model with higher 
magnetic field strength, all simulations produce contrast values 
for the undegraded images that lie within their mutual varia- 
tion ranges (see right column in the figure). After degradation 
(middle column), however, a systematic shift can be seen. The 
C0 5 BOLD models are stacked according to the size of their spa- 
tial extent: low contrasts for model by Steffen to higher contrasts 
for the smallest model Cm. The systematic difference can be re- 
moved by consistently going towards a large FOV of the same 
size. Making a realistic extension of the FOV is only practically 
feasible for disc-centre images. For regions towards the limb, 
however, the extension of the FOV would require to vary the ob- 
serving angle across the image. As there is no straightforward 
procedure to do this, we refrain from enlarging the PSF for the 
calculation of the centre-to-limb variation of the contrast. 

The uncertainty in the PSF wing contribution (as seen from 
the choice of the FOV) is counteracted by uncertainties in the 
PSF parameters. In the following, we illustrate the effect for the 
blue continuum images. So far we used the best-fit PSF as de- 
termined in Paper I. The PSF for the blue BFI channel is deter- 
mined by the Voigt parameters y = 0'.'004 and cr = 0"008 (see 
Sect. 12.21 ). Application of this PSF produces mean contrasts of 
13.0 % and 1 1.3 % for the small and the large FOV, respectively. 
We now consider the error margins for the non-ideal Voigt con- 
tribution to the PSF as given in Table 4 in Paper I. First we re- 
main with y = 0'.'004 and shift cr according to the uncertainties 
of the linear regression stated in Table 4 in Paper I. The resulting 
mean contrast range from 12.9% to 13.4% for the small FOV 
to 1 1.2 % and 1 1.8 % for the large FOV. The spread in possible 
contrast values becomes larger if also the uncertainty in the pa- 
rameter y is considered. Within the specified error ranges, we 
can produce contrasts up to 13.2% or even 13.6%. The intrin- 
sic uncertainty in the PSF apparently allows for the contrasts 



of the degraded synthetic maps in the range of ~ 11.2% to 
~ 13.6 %. The observational result of 12.8 % is well included in 
this range. In return, the uncertainty in the non-ideal PSF contri- 
bution would also lead to corresponding error margins for the de- 
convolved images so that also the original synthetic images and 
the deconvolved observations agree in terms of contrast within 
their mutual error margins. 

Limb observations. The solar limb is used to re-determine the 
coordinates of the FOV whenever it can be seen in an image. The 
derived positions sometimes differ from the Hinode coordinates 
in the image file headers by of the order of 1 " or less. At the limb, 
a small pointing error can translate into a significant uncertainty 
in fi (e.g.: A/j. « 0.05 for 1"). The unavoidable smearing due to 
the limited spatial resolution of an instrument limits the accu- 
racy of the limb determination. The uncertainty can be given as 
A/i < 0.03 for the BFI continuum bands, when assuming one 
FWHM of the PSF core as maximum offset. Effects like thermal 
deformation of the instrument structure and fluctuations in the 
flight characterist i cs of th e spacecraft induce an additional jitter- 
ing. iKosugi et al.l {2007) state that the image stabilization sys- 
tem of SOT achieves 0"007 (3 cr level) and effe ctively restricts 
the rms image displacem ent to less than 0"03 dTsuneta & et al.l 
120071: iKosugi et al. 2007). The influence of the spacecraft on the 
spatial resolution should thus be rather insignificant. We never- 
theless explicitly state that data for ju < 0.1, corresponding to 
distances less than ~ 5 " from the limb, has to be regarded with 
great care in view of the potential error margins. 

Power excess and JPEG compression. In Sect. 12.31 we men- 
tioned a small power contribution at spatial scales between the 
granulation peak an d the spatial resolution of SOT (see Fig.|9]l. 
Wog er et al.l (|2008) describe a similar power excess between 
SOT and DST data (see their Fig. 5). It is clearly in excess com- 
pared to synthetic images and also is usually not seen in compa- 
rable SST data. We therefore argue that the power excess is an 
artefact. As the enhanced power peak occurs at frequencies that 
correspond to roughly 8 pixels and less (see Fig. |9), it could in 
principle be related to JPEG compression. Essentially all SOT 
images are compressed before downlinking to the ground in or- 
der to increase the telemetry throughput. The used lossy JPEG 
image compression uses a 8 x 8 pixel blocking. Consequently 
artefacts can be expected on pixel scales of 8 and harmonics 
thereof. We therefore compressed and decompressed synthetic 
images with routines kindly provided by Katsukawa & Tsuneta 
(priv. comm.). The power spectral density of the resulting pro- 
cessed images do not significantly differ from those of the orig- 
inal images, disqualifying the JPEG compression as the source. 
Whatever the cause of the artefact, the effect on the intensity con- 
trast is negligible as the affected scales do not contribute much 
to the overall contrast. 

8. Conclusion 

State-of-the-art radiation (magneto-)hydrodynamic simulations 
can reproduce many properties of the observed continuum inten- 
sity distribution if the effects of instrumental image degradation 
are properly taken into account. The latter can have dramatic 
effects. Although SOT is an excellent instrument, which works 
close to the diffraction limit, the non-ideal component to the PSF 
has a significant effect on the intensity distribution. The removal 
of the degradation by PSF-deconvolution more than doubles the 
intensity contrast. In contrast to earlier studies, we derive final 
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contrast values at blue, green, and red continuum wavelengths 
that are slightly higher than for the original synthetic images, 
although they still agree within the mutual variation ranges. A 
detailed realistic treatment of instrumental straylight is found to 
be of crucial importance for a quantitative comparison between 
simulations and observations. A constant straylight contribution 
alone is not sufficient and misleading. 

Owing to the many sources of uncertainty, the intensity con- 
trast as a single number is a poor mean of comparing observa- 
tions and simulations. This is already problematic for compar- 
isons of observations that are taken with different instruments or 
under different conditions. The only practicable way is to prop- 
erly remove the effects of image degradation prior to comparison 
as it is demonstrated here for SOT. This procedure is unfortu- 
natly not always possible as the properties of the PSF are often 
poorly known. However, even in such rather fortunate cases as 
for SOT, one cannot expect an exact match between observa- 
tions and simulations. Already the uncertainties in the non-ideal 
PSF contributions, result in a error margins of the order of 1 or 
possibly even 2 percent points in the blue continuum. 

Comparisons should therefore be based on more informative 
properties of the intensity distribution. Here we find the intensity 
histograms and the power spectral density from observation and 
simulation to agree very well. It is obvious that all spatial scales 
with significant contributions to the contrast are resolved with 
the 50-cm aperture of SOT. Therefore we do not expect that ob- 
servations with larger aperture telescopes would yield different 
contrast values for the quiet Sun. We finally conclude that the 
traditionally perceived conflict between observations and simu- 
lations can now be dismissed. 
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